Electronic transport in quasi-one-dimensional arrays of gold nanocrystals 



in 

O 

o 

(N 
oo 

(N 



Klara Elteto 1 , Xiao-Min Lin 2 , and Heinrich M. Jaeger 1 
1 James Franck Institute and Department of Physics, University of Chicago, Chicago, IL 60637 
2 Materials Science Division, Chemistry Division and Center for Nanoscale Materials, 
Argonne National Laboratory, Argonne, IL 60439 
(Dated: February 2, 2008) 

We report on the fabrication and current-voltage (IV) characteristics of very narrow, strip-like 
arrays of metal nanoparticles. The arrays were formed from gold nanocrystals self-assembled between 
in-plane electrodes. Local cross-linking of the ligands by exposure to a focused electron beam and 
subsequent removal of the unexposed regions produced arrays as narrow as four particles wide and 
sixty particles long, with high degree of structural ordering. Remarkably, even for such quasi-one- 
dimensional strips, we find nonlinear, power-law IV characteristics similar to that of much wider 
two-dimensional (2D) arrays. However, in contrast to the robust behavior of 2D arrays, the shape 
of the IV characteristics is much more sensitive to temperature changes and temperature cycling. 
Furthermore, at low temperatures we observe pronounced two-level current fluctuations, indicative 
of discrete rearrangements in the current paths. We associate this behavior with the inherent high 
sensitivity of single electron tunneling to the polarization caused by the quenched offset charges in 
the underlying substrate. 

PACS numbers: 73.23.-b, 73.22.-f, 73.23.Hk 
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I. INTRODUCTION 

Metal or semiconductor nanoparticles and complex 
structures built from them through self-assemblyiS*^ 
provide model systems for further investigation of physics 
on the mesoscopic scale, where quantum confinement and 
Coulomb charging set the relevant energy scales^ Elec- 
tronic properties of single particles^ as well as two- and 
three-dimensional superlattices of nanocrystalsi&Si have 
been extensively studied over the last decade. Here we 
address what happens to electronic transport as the one- 
dimensional limit is approached. 

Electronic transport through an array of small metal 
particles separated by nanoscale gaps is determined by 
the interplay between single electron charging of an indi- 
vidual particle and tunneling between adjacent particles. 
In the presence of charge disorder due to quenched impu- 
rities in the insulating substrate under the array^ this 
interplay leads to highly non-Ohmic current-voltage TV 
characteristics. Theory and simulations predict that, at 
sufficiently low temperatures, no current flows below a 
voltage threshold, Vt, while above it the current follows 
a power-law, / ~ (V — Vtp , with £ = 1 in. ID and 
between 5/3 and 2 in 2D. 10 i 1:L i 12 

Experimentally, true ID chains are difficult to achieve 
and, so far, results have been obtained only on quasi- 
1D structures with significant amount of disorder in 
the particle arrangements. Narrow chains of carbon 
nanoparticles^ showed sample-dependent scaling expo- 
nents between 1 and 2.3. Electron-beam written lOOnm 
wide multilayers of AU55 particles^ exhibited £ = 1.6, 
and experiments on strips of 1.8nm gold clusters^* gave 
£ = 3. This spread in exponents is in sharp contrast 
with the situation for well-ordered, 2D Au monolayer 
nanocrystal arrays. Experiments on highly ordered 2D 
nanocrystal arrays found C = 2.25 ±0.1, a value that is 



temperature-independent and highly reproducible from 
sample to sampleAi& 

These results highlight two key issues. First, to what 
extent is the spread in the measured scaling exponents 
caused by the charge disorder in (quasi-) ID structures 
rather than to structural disorder and fabrication details; 
second, how narrow does an array have to be in order 
for its current-voltage characteristics to cross over from 
a nonlinear, 2D behavior to the a linear response with 
£ = 1? To address these issues, we performed system- 
atic measurements on quasi- ID strips fabricated from 
well-ordered monolayer superlattices. These structures 
also allow us to observe explicitly the effects of quenched 
charges in the substrate through measuring the current 
fluctuation at fixed temperature and the changes in the 
current-voltage characteristics after temperature-cycling 
the array. 



II. SAMPLES 

Gold nanocrystals of 5.5nm diameter with 5-7% dis- 
persity were synthesized as described in Ref. ^3 through 
a digestive ripening method and suspended in toluene. 
A droplet of this colloidal suspension was then de- 
posited onto Si3N4-coated Si substrates with prefabri- 
cated, 15nm thick in-plane chromium electrodes (2/xm 
wide with 500nm gap, i.e., about 60 particles across). By 
controlling the particle concentration, solvent evapora- 
tion rate and concentration of excess alkanethiol ligands 
in the solution, this drop drying procedure was used to 
form highly ordered nanocrystal monolayers^ Alterna- 
tively, a water droplet was first deposited onto the sub- 
strate and the colloidal suspension was allowed to creep 
up at the water-air interface, forming a thin layer which 
then draped itself over the substrate as the water slowly 
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evaporatedJ-Si Both methods allow nanocrystals to self- 
assemble into a compact 2D monolayer at the liquid-air 
interface before the solvent dewetting occurs. 

After solvent evaporation, the quasi-lD array was fab- 
ricated by exposing the corresponding monolayer regions 
to a finely focused electron beam in a scanning electron 
microscope (SEM)i 20 i 21 i 22 Line doses of 8 — 25nC/cm or 
area doses of 10mC/cm 2 were used to e-beam- write lines 
or rectangles at 30kV. These exposure parameters effec- 
tively cross-link the interstitial ligands but do not affect 
the nanoparticles' shape and relative position^ The ex- 
posed regions become resilient to heated toluene, which 
can be used to wash away unexposed particles. Overall, 
this method works like a positive e-beam resist in that the 
final structures correspond to the exposed areas. Back- 
etched SiaN4 window areas in the substrate allowed for 
inspection by transmission electron microscopy (TEM). 
In our experiments, all the samples were inspected by 
TEM after completing the electronic measurements. 

Figure ^ shows TEM images of structures fabricated 
by this method. The shape of the final patterns depends 
only on the ability to control the position and movement 
of the electron beam, while the ultimate resolution ap- 
pears to be limited primarily by the e-beam spot size and 
the degree of ordering in the original monolayer. Panel 
a gives the details of a narrow strip fabricated through 
this method. Panel b shows one of the measured arrays 
with the Cr electrodes visible. All samples were struc- 
turally ordered in the sense that they retained the close- 
packed particle arrangement of the monolayer, showed 
few defects, exhibited a uniform inter-particle spacing, 
and had well-defined edges with a roughness of 1-2 parti- 
cles (in many cases unavoidable because of the mismatch 
between the particle lattice orientation and the direction 
of e-beam writing). 

Six arrays were beam-written as single lines and had 
average widths of 30nm, or four particles across. Two 
arrays were patterned as slightly wider rectangles, on av- 
erage seven particles across. Isolated, localized patches 
where particles formed a second layer made up less than 
15% of the area in each of the eight arrays. 

All samples were mounted in a shielded cryostat and 
connected with low-noise coaxial cables to a Keithlcy 
6340 sourcemeter. Temperature was controlled from 8K 
to 120K, with stability better than 0.1%. We measured 
the current while ramping the bias voltage between -20V 
and +20V at rates of 50mV/sec or less. The rms current 
noise of the set-up was about lOfA. 

As a control study, we measured full 2D monolayers 
before and after e-beam writing and washing. We found 
that the voltage threshold at a given temperature, and 
the scaling exponent did not change by more than the 
fitting uncertainty. This shows that the alkanethiol lig- 
ands act as mechanical spacers and, cross-linked or not, 
do not otherwise influence the properties of the tunneling 
barrier between the neighboring particles. 
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FIG. 1: Transmission electron micrographs of structures pat- 
terned from an extended nanoparticle monolayer using elec- 
tron beam writing, a) Large magnification image of a nar- 
row strip, showing the integrity of the particle arrangement 
after e-beam exposure, b) 30nm wide, quasi-lD array be- 
tween planar Cr electrodes. Gaps between the particles are 
1.3nm on average, c) Beam- written array formed using a col- 
loid containing excess dodecanethiol ligand. The array folds 
over, as shown here, rather than break apart during aggres- 
sive agitation in heated toluene, indicating e-beam-induced 
cross-linking of the ligands. d) Arbitrary patterns with both 
straight and round sections are possible. 



III. RESULTS AND DISCUSSION 

Typical current-voltage curves at different tempera- 
tures are shown in Fig. Data from all measured 
arrays follows power-law scaling, / ~ (V— Vt)^, with volt- 
age thresholds, Vt, in the several- volt regime and scaling 
exponents, £, around two. Figure |2Jd shows the extent 
of the power-law scaling. For each temperature, we de- 
termined the voltage threshold and the scaling exponent 
from best straight-line fits on a log-log plot. On closer 
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FIG. 2: Current-voltage (IV) characteristics of a 4-particle- 
wide array, a) The shape of the IV curves changes with tem- 
perature, unlike what is observed for wider, 2D arrays, b) 
Log-log plot of the IV curve for the same array as in (a) at 
20K with best-fit voltage threshold V*(20K) = 2.5V. The line 
corresponds to a scaling exponent £ = 1.95. c) Temperature- 
dependent values of £ for all eight samples measured. The 
hatched area gives the standard deviation in the observed 
values for £ in 2D arrays. 



inspection, Fig. [2^ also shows that the shape of the IV 
traces changes with temperature. This is in contrast to 
what has been observed for wider, 2D arrays, where the 
scaling exponent is temperature-independent and the IV 
curves can be collapsed simply by a translation along the 
voltage axis^ In Fig. we plot both the extent of the 
temperature-induced variation in £ for a given sample 
and the variation of C for different samples. For compar- 
ison, we also show the very small variation of the scaling 
exponent obtained from a series of 2D arrays fabricated 
with the same Au nanoparticles^ 

For Vt, on the other hand, we find behavior similar to 
what has been previously observed in 2D arraysi&2£ and 
in ID chains^ namely an approximately linear decrease 
with increasing temperature. At zero temperature, the 
voltage threshold can be expressed as V t (0) = aNVo, 
where N is the number of particles along the length 
of the array and eVb corresponds to the typical electro- 
static energy cost associated with single electron tunnel- 
ing between neighboring particles in the array^ Because 
quenched charge disorder leads to a distribution of lo- 



cal energy costs, tunneling will occur along those current 
paths that minimize penalties and circumvent sites with 
large local charging energies. The extent to which this 
happens is measured by the prefactor a which depends 
on the array geometry and dimensionality. 

For 2D close-packed arrays, characterized by a large 
coordination number and thus a multitude of possible 
detours, theory predicts a = 0.226, while a = 1/2 for 
ID chains>i2ii& To find a experimentally, we obtained N 
directly from TEM images and established Vt(0) by ex- 
trapolating the temperature-dependent threshold to zero 
temperature. Using the procedure described in Ref. 0, 
Vo was calculated from the average radius and center- 
to-center spacing of the particles in each array, as de- 
termined by TEM. Our quasi- ID strips exhibit values 
between the ID and 2D theoretical prediction, giving 
a = 0.34 ± 0.09 for the eight samples measured. 

The observed values for a and £ demonstrate that, 
even at average widths of only four particles, transport 
along the arrays has not yet reached the ID limit. The 
nonlinear shape of IF curves shows that the current paths 
are still able to meander and branch significantly, which 
is the characteristic behavior of 2D transport. 

Nevertheless, quasi- ID arrays allow for only a very lim- 
ited number of parallel paths. We therefore expect that 
these limited number of current-carrying paths will be in- 
herently more susceptible to changes in the configuration 
of quenched offset charges in the underlying substrate. 

This is borne out by measurements of the fluctuations 
in the transport current at low temperatures under a 
fixed bias voltage. We observe clear evidence for ran- 
dom switching between a few well-defined states of the 
system. The I(t) trace at 8K in Fig. [21} shows that 
the current fluctuates on the order of several minutes be- 
tween two steady states, which also are easily identified 
from the histogram in Fig. |3p. The magnitude of the 
current change in this figure, about 12%, corresponds to 
one out of four parallel conduction channels in the array 
changing its transmission by 50%. 

As the temperature increases between 8K and 60K, the 
frequency of this telegraph noise increases and, above 
60K, distinct current states can no longer be observed 
(Fig. 13;). Qualitatively, this change in behavior is ex- 
pected if the number of active conduction channels in 
the array depends on the configuration of quenched off- 
set charges in the substrate. With increasing thermal 
energy, these offset charges can switch more rapidly be- 
tween different trapping states, inducing the opening and 
closing of conduction paths. Similar stochastic switching 
of current was also observed in much wider 2D arrays 
(with width of 2/xm). However, these 2D arrays did not 
exhibit bimodal switching even at 8K. 

Further evidence that rearrangements of trapped offset 
charges determine the behavior of the mobile electrons 
comes from temperature-cycling experiments. In these 
experiments, we measure the current at 20K, warm up 
the array to 200K for at least an hour and, and upon 
cooling back to 20K, measure the current again (Fig. 
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FIG. 3: Current fluctuations at fixed bias voltage, a) Per- 
centage change in current as a function of time for an array 
at 8K and -12V bias. At the large jumps, the current value 
changes by as much as lpA between successive data-points 
(1 second), b) Histogram of the data in (a), c) Percentage 
change in current as a function of time for the same array at 
100K and -12V bias. 



Et). Such temperature cycles are done several times for 
the same sample. We find that, in response, the thresh- 
old voltage and scaling exponent change randomly (Fig. 
0|d,c), varying significantly more than the fitting uncer- 
tainty given by the error bars. In fact, the changes in- 
duced by temperature cycling are of similar magnitude 
as the sample-to-sample variations (Fig. |2t). This vari- 
ability contrasts with the robust behavior we find in 2D 
arrays that are a few hundred particles wide. For 2D ar- 
rays, even after temperature cycling, £ stays within the 
range indicated in Fig. |2J;. 
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FIG. 4: The effect of temperature cycling, a) Current as a 
function of voltage for a 4-particle-wide array at 20K after 
several temperature cycles to 200K. b) Voltage threshold, Vt, 
of the array obtained from power-law fitting of IV curves, c) 
Scaling exponent, £ , of the IV curves. The error bars indicate 
the fitting uncertainty. 



the quasi-lD nature of the samples. The fact that there 
are fewer conduction paths in narrow strips as compared 
to the extended 2D arrays implies a high sensitivity to 
the change of the quenched charge disorder "landscape" . 
This charge reconfiguration is time dependent and in- 
duced by thermal activation. 

In a 2D array by contrast, changes in the charge 
configuration and tunneling paths average out, giving 
a time- and temperature-independent scaling exponent 
that varies little from array to array. Thus, while 
quenched charge disorder gives rise to the nonlinear 
current-voltage characteristics in the first place, quasi- 
lD structures at any given time sample only a limited 
range of disorder configurations. This has a significant 
implication for device applications in which a high degree 
of reproducibility from sample to sample is desired. In 
those cases, it will be necessary to use wider nanocrystal 
arrays in which there is sufficient spatial averaging over 
different trapped charge disorder configurations. 
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IV. CONCLUSIONS 

In summary, our results demonstrate the important 
role played by quenched or trapped offset charges in de- 
termining the arrays' overall transport behavior. Being 
able to control the structural integrity of the arrays down 
to widths of 4 particles, we conclude that the large vari- 
ability in threshold voltages and scaling exponents seen in 
our experiments as well as by othersi^i^i^ is inherent to 
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